Common trends in the critical behavior of the Ising and directed walk models 
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We consider layered two-dimensional Ising and directed walk models and show that the two problems 
are inherently related. The information about the zero-field thermodynamical properties of the 
Ising model is contained into the transfer matrix of the directed walk. For several hierarchical and 
aperiodic distributions of the couplings, critical exponents for the two problems are obtained exactly 
through rcnormalization. 
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The Ising model (IM) and the directed walk (DW) are 
among the most studied problems in lattice statistics. 
The IM is a standard model for magnetic or liquid-gas 
phase transitions whereas the DW can be used to describe 
linear fluctuating objects such as directed polymers, flux 
lines or interfaces in two-dimensional systems. 

The IM is exactly solvable in two dimensions 0] and 
the solution can be generalized for layered systems with 
different types of distributions for the intcrlaycr couplings 
such as periodic g] , quasi-periodic || , aperiodic and 
random || . The DW is probably the simplest non-trivial 
problem in statistical mechanics for which exact results 
can be obtained on homogeneous J?| , inhomogeneous ||] 
and random [^| lattices. 

In this Letter, we present a hitherto unnoticed con- 
nection between the IM and the DW in two dimensions. 
Both problems are considered on layered lattices, such 
that the walk is directed along the translationally invari- 
ant direction. We show that the complete solution of the 
DW, i.e. the diagonalization of its transfer matrix (TM), 
provides all the necessary information to obtain the zero- 
field thermodynamical properties and correlation func- 
tions of the IM. The DW approach, which is simpler, 
is then used to perform an exact renormalization-group 
(RG) study of the TM eigenvalue problem for self-similar 
distributions of the couplings. The critical properties of 
the IM and DW are governed by two different fixed points 
of the same RG-transformation. 

Let us first present the hidden relation between the 
two problems. We consider a layered IM in the extreme 
anisotropic limit flj0| . The transfer matrix going in the 
direction parallel to the layers is exp(— tH), where r is 
the lattice spacing in the Euclidian time direction, and 
TL the Hamiltonian of a quantum Ising chain: 
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The (jjS are Pauli spin operators, the transverse field hk 
h plays the role of the temperature and the couplings J k 
are non-periodic. 

Following Lieb et al and Pfeuty |Q H can be 
transformed into a free-fermion model 
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in terms of the fermion creation and annihilation opera- 
tors 77J, rj q . The fermion excitations A ? are non- negative 
and satisfy the set of equations 

Aq^q(k) = -h k $ q (k) - J k <S> q {k + 1) 

A 9 $ 9 (fc) = - J k -i^ q (k - 1) - h k ^ q {k) (3) 

with the boundary conditions Jq = Jl = 0. The & q s and 
^ q s, which are related to the coefficients of a canonical 
transformation, are normalized. They enter into the ex- 
pressions of correlation functions and thermodynamical 
quantities 

Usually one proceeds by eliminating either *& q or *& g 
in ([}]) and the excitations are deduced from the solution 
of an eigenvalue problem. This last step can be avoided 
by introducing a 2L-dimensional vector V g with compo- 
nents V q (2k-l) = -$q(k), V q (2k) = * g (fe) and noticing 
that the relations in Eq. (||) correspond to the eigenvalue 
problem for the matrix 
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which can be interpreted as the TM of a DW problem on 
two interpenetrating, diagonally layered square lattices. 
The walker makes steps with weights h k and J k between 
first-neighbour sites on one of the two lattices. 

Changing <S> q into — 3> 9 in V g , the eigenvector corre- 
sponding to — A q is obtained. Thus all the information 
about the DW and the IM is contained into that part 
of the spectrum with A q > 0. Later on we shall restrict 
ourselves to this sector. 

Let us now consider the correlation lengths in the di- 
rection parallel to the layers for both problems. For the 
DW it can be expressed as a function of the two leading 
eigenvalues of the TM with: 
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FIG. 1. Matrix-elements Tjj+i as a function of j for differ- 
ent non-periodic sequences: a) hierarchical b) period-doubling 
c) three-folding d) paper-folding. Components of the eigen- 
vector to be decimated out in the RG transformation are de- 
noted by crosses. The heights of solid vertical bars indicate 
the strength of the couplings, the grey bars stand for the field. 

Approaching one of the two critical points, the cor- 
relation length of the problem is diverging and the cor- 
responding part in the TM spectrum displays a scaling 
behavior. Let us consider a finite system with transverse 
size L ;§> 1 and denote by AA^ either A^ — A^_i for the 
DW or Ai itself for the IM with i < L. When lengths 
are rescaled by a factor of b > 1, i.e. with L' = L/b, the 
gaps are assumed to behave as 
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where y\ is the gap exponent which is generally different 
at different parts of the spectrum. This leads to the finite 
size behavior 



AAj(L) ~ L~ 



(8) 



thus from Eqs. (|5|) and (0) the longitudinal correlation 
lengths are £ii ~ L VA . Since £j_ ~ L, the anisotropy 
exponent z, such as ^ ~ is given by z = y\. For the 
DW one is interested in the transverse fluctuations which 



are characterised by the wandering exponent w through 
£_L ~ fjj", thus w = yl 1 . 

The scaling properties of the spectrum of T are conve- 
niently studied using RG techniques. We consider differ- 
ent self-similar lattices for which exact RG transforma- 
tions can be worked out so that we obtain exact results 
about the critical properties of both the IM and the DW. 
In the following the transverse field is assumed to be con- 
stant and equal to h. 

Hierarchical sequence We start with a hierarchical 
lattice in which the couplings follow the Huberman- 
Kerszberg sequence |Q , 



J k = R n J, fc = 2 n (2m+l). 



0,1, 
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with < R < 1. The eigenvalue problem for T corre- 
sponds to the second order difference equations 

T^VU - 1) - A.V (J) + T jd+1 V(j + 1) = , (10) 

where j = 1, • • • , 2L. To construct an exact recursion we 
eliminate from these equations components of the form 
V{41 + 2), V(Al + 3) which are connected to a J coupling 
(indicated by crosses in Fig. la). After such a decimation 
the triplet (h, J, h) is replaced by a renormalized field 
h! and keeping R unchanged, the remaining couplings 
become J' — RJ due to the hierarchical structure of 
the sequence. Thus the renormalized equations keep the 
original form with A changed into A'. Introducing the 
reduced variables A = J/h and A = A/h one arrives at 
the two-parameters recursion: 



A' = y (A 2 — A 2 — 1) 
A 



A' = R (A 2 - A 2 ) . (11) 



The RG-transformation has two non-trivial fixed points, 
governing the scaling of the eigenstates at the top of the 
spectrum (DW) and at A = (IM) , respectively. The line 
A = 0, corresponding to the IM situation, is invariant 
under the RG transformation in Eq. (|ll]). Along this 
line, starting with a ferromagnetic model with A > 0, 
after one recursion step the system is transformed into 
an antiferromagnetic model with A < 0. The critical IM 
with A c = l/R is transformed into the IM fixed point, 
which is situated at A* = —l/R. At the IM fixed-point 
the leading eigenvalue of the transformation is ei = R + 
1 /R and the anisotropy exponent of the hierarchical IM 
is given by: 



z = VA = 
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Thus scaling in the hierarchical IM close to the critical 
point is essentially anisotropic. 

Scaling of the eigenstates at the top of the spectrum is 
governed by the DW fixed-point situated at 
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and the leading eigenvalue is given by: 



ei 



1/2 



Thus the wandering exponent of the walk is: 
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In the homogeneous model, with R = 1, the DW fixed 
point is shifted to infinity since A* — » oo, A* — > oo and 
along the separatrix A/A — ► 1. To evaluate the scaling 
behavior we introduce new variables: k = 1/A and A = 
(A/A) 2 — 1, in terms of which the fixed point is given by 
k* = and A* = 0. Then the separatrix is a straight 
line: A(k) = a*n, with a* = 2, and according to Eq. 
one point of the (k, A) plane with A = ok will transform 
into («' = k/cl, A' = (1 - 2/a 2 )A). Thus a' = a 2 - 2 and 
the leading eigenvalue of the transformation is e% = 4, 
consequently w = 1/va = 1/2, in agreement with known 
results 0. We note that w(R) is discontinuous at R = 1, 
since from Eq. @ lim fl _>i_ w(i?) < 1/2. 

Period- doubling sequence In our next example, the 
couplings Jfc are generated according to the period- 
doupling sequence |l4| which follows from the substitu- 
tion A — > AB and B — > AA. Here and in the following, 
the couplings are parametrized as J a = J and Jb = RJ ■ 

In an exact RG transformation, six sites out of eight 
have to be decimated, as indicated on Fig. lb. Asso- 
ciating new couplings h' with the decimated blocks one 
obtains a recursion in terms of h! and A' while A and 
R\ thus the ratio i?, remain unchanged. In terms of the 
reduced parameters the RG-transformation reads as 
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with c = A 2 (-A 2 + 1 + A 2 ) 2 - i? 2 A 2 (A 2 - A 2 ) 2 and 
d = (K 2 - l) 2 - A 2 A 2 (1 + R 2 ) + A 2 (l + i? 2 A 2 ). The 
IM-fixed point of the transformation is at A* = and 
A* = -i? -1 / 3 , with the leading eigenvalue ei = (i? 1 / 3 + 
Since the rescaling factor of the transformation 
is 6 = 4 we obtain 
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for the anisotropy exponent of the period-doubling IM. 

The top of the spectrum, corresponding to the DW 
problem, scales to a fixed point with A — * oo, A — * oo, 
but A/A — ► R. In terms of the variables K = 1/A 
and A = (A/A) 2 - R 2 the fixed-point is at n* = 
and A* =0, while the separatrix, close to the fixed 
point, is of the form A(k) = a*n 2 + 0(k 4 ), with a* = 
(V2R - 2R 2 )/{1 - R 2 ). Then, according to Eq. ©, a 
point of the (k, A)-plane with A = an 2 will transform 



to (k' ~ k 2 , A' ~ A 1 / 2 ). This type of scaling behavior 
is compatible with an essential singularity in the gaps at 
the top of the spectrum, 



AAi - exp(-CL' T ) 



(18) 



with a = 1/2, since the rescaling factor is b = 4. Thus the 
parallel correlation length of the DW is given by Cy >W ~ 
exp(CL 1 / 2 ) and the transverse fluctuations of the walk 
grow anomalously, on a logarithmic scale: 



[x(t)-x(o)y 



1/2 



ln 2 W 



(19) 



Here X(t) denotes the position of the walker at time t. 
We note that the same asymptotic behavior is found in 
the Sinai model |l5| of a one-dimensional random walk 
in a random environment. 

Three-folding sequence The three-folding sequence is 
generated by the substitutions A — * ABA, B — * ABB 
|]l6| . In the RG transformation - as indicated on Fig. lc - 
blocks of four sites are decimated out. Due to the asym- 
metric nature of the blocks, after one RG step the trans- 
fer matrix becomes asymmetric, too: Tj j+i/Tj+i j = s 
for j even, while Tjj + i/Tj + ij = s^ 1 for j odd. 

The recursion relations in this case are more conve- 
niently expressed using the variables A = A/J, /i = h/J 
and s, while R remains unchanged: 
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with c = M 2 (A 2 - y 2 - R 2 ), d = /i 2 (A 2 - n 2 - 1) and e = 
(A 2 - 1)(A 2 - R 2 ) - ^ 2 A 2 . We note that the asymmetry 
parameter s, such that s' = s(c—e)(d—e), does not enter 
into the recursions for A and fi. 

At the IM fixed point (A* =0, n* = R 1/2 ) the leading 
eigenvalue of the RG transformation is t\ = [(2 + R)(2 + 
i? -1 )] 1 / 2 , thus the anisotropy exponent is given by 



ln(2 + R)(2 + R~ 
2 In 3 



(21) 



The DW fixed point is again at infinity: A* = oo, 
/i* = oo, with A* / fi* = 1. The scaling behavior at this 
fixed point is similar to that in the period-doubling case. 
The eigenvalues at the top of the spectrum show an es- 
sential singularity like in Eq. ( |l8| ) with a = 1/2 and the 
transverse fluctuations grow on a logarithmic scale as in 
Eq. ©. 

Paper-folding sequence Finally, we consider the pa- 
per-folding sequence p6| which is generated by the two- 
letter substitutions AA -> AABA, AB -> AABB, 
BA -> ABBA and BB -> ABBB. In the RG trans- 
formation, decimating out blocks of two sites (Fig. Id), 
alternating field variables hi and hi are generated for 
odd and even lattice sites, respectively. Furthermore, the 
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transfer-matrix becomes asymmetric and the asymmetry 
parameters are different for odd and even elements. As a 
consequence, the exact RG transformation contains alto- 
gether six parameters. Here we just present the scaling 
behavior at the two non-trivial fixed points, details of the 
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calculation will be presented elsewhere |17|. 

At the IM fixed point, the anisotropy exponent is con- 
tinuously varying and given by: 



+ R){1 + R- 1 ) 
In 4 
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At the DW fixed point the scaling is again of the streched 
exponential form with a leading behaviour for transverse 
fluctuations given by Eq. (|u}). 

Let us now turn to a discussion of the critical behavior 
we have obtained for the IM and the DW. All the ape- 
riodic IMs we considered are strongly anisotropic with 
a continuously varying anisotropy exponent. In the ex- 
tended parameter space there is a line of fixed points 
parametrized by the coupling ratio R. The critical be- 
havior of the DWs is also found to be anomalous: the 
lines of fixed points of the non-periodic systems are dis- 
connected from the fixed point of the homogeneous ones. 
For the hierarchical model the wandering exponent is dis- 
continuous at R = 1 whereas for the other sequences 
the transverse fluctuations grow on the same logarithmic 
scale. 

The difference between the IM and the DW on 
the same lattice can be understood using a relevance- 
irrelevance criterion |Q| which is a counterpart for ape- 
riodic systems of the Harris criterion Jl8| for random 
ones. The cross-over exponent associated with a layered 
non-periodic perturbation is (f> = 1 + — 1) where 
v is the exponent of the correlation length, 'perpendicu- 
lar to the layers, for the unperturbed system and f2 is 
a wandering exponent |19| which characterizes the fluc- 
tuations in the couplings Jk around their average J as 

All the non-periodic sequences we considered have 
n = 0. For the IM with v = 1 the cross-over expo- 
nent vanishes. Thus the perturbation is marginal, which 
explains the continuous variation of the anisotropy expo- 
nent. On the other hand, for the anisotropic DW with 
2/_L = 1/2, the cross-over exponent is 4> — 1/2 and the 
perturbation is relevant. This is again in agreement with 
our results. 
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